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Abstract. Integrable perturbations of the nonholonomic Suslov, Veselova, Chaplygin and 
Heisenberg problems are discussed in the framework of the classical Bertrand-Darboux 
method. We study the relations between the Bertrand-Darboux type equations, well stu¬ 
died in the holonomic case, with their nonholonomic counterparts and apply the results to 
the construction of nonholonomic integrable potentials from the known potentials in the 
holonomic case. 
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1 Introduction 


In classical mechanics, the Euler-Poisson equations 

( 1 - 

describe the rotation of a rigid body with a fixed point using a rotating reference frame with its 
axes fixed in the body and parallel to the body’s principal axes of inertia. Here oj = (uq, W 2 , o; 3 ) is 
the angular velocity vector of the body, I = diag(/i, I2,13) is a tensor of inertia, 7 = ( 71 , 72,73) is 
a unit Poisson vector and V ( 7 ) is a potential field. All the vectors are expressed in the so-called 
body frame and x x y means the cross product of two vectors in three-dimensional Euclidean 
space. 

Let us impose a nonholonomic constraint on the angular velocity 


Id; = Iw x oj + 7 x 


d'y 


7 = 7 x 0 ; 


/ = ( 0 ;, a) = 0 or f = (uj, 7 ) = 0 , 


where a is a fixed unit vector in the rotating frame for the Suslov problem [33], 7 is a fixed unit 
vector in the stationary frame for the Veselova problem [40] and ( 2 , y) means the scalar product 
of two vectors. In this case the Euler-Poisson equations (1.1) are replaced by equations 


Iw = Ixw + 7 X 


dv(rt) 

d'y 


+ An, 


7 = 7 x 0 ;, 


n = a, 7 , 


( 1 . 2 ) 


where A is a Lagrange multiplier which has to be found from the condition / = 0. 

*This paper is a contribution to the Special Issue on Analytical Mechanics and Differential Geometry in honour 
of Sergio Benenti. The full collection is available at http://www.emis.de/journals/SIGMA/Benenti.html 
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Both systems of differential equations (1.1) and (1.2) are geometrically interpreted in terms 
of the vector field X 

±i = Xi{x i,...,x 6 ) (1.3) 

in a six-dimensional manifold M with coordinates x = (a;, 7 ). The classical Euler-Jacobi theorem 
says that the vector field X (1.3) on a six-dimensional manifold M is integrable by quadratures 
if it has an invariant volume form (invariant measure) and four functionally independent first 
integrals [28, 37]. 

Equations (1.1) and (1.2) preserve the norm of the unit Poisson vector 7 

Ci = ( 7 ) 7 ) = ll + 72 + 73 = !> 
and mechanical energy 

Hi = ^(Jw,w) + ^(7). 

An additive perturbations (1.2) of the Euler-Poisson vector field (1.1) change the standard 
invariant volume form and the second geometric first integral 

C 2 = (oj, I7) -7 C 2 = (u , n), 

see details in [15]. Thus, according to the Euler-Jacobi theorem equations (1.1) and (1.2) are 
integrable by quadratures if there is one more independent first integral H 2 . 

There are several methods to uncover integrals of motion. A one of the simplest method 
considers function H 2 polynomial in the velocities, with the coefficients being arbitrary functions 
of the coordinates, and requires that it is conserved in time 

H 2 = 0 . (1.4) 

This condition yields a system of coupled partial differential equations on the coefficients, which 
are some of the most well-studied first-order PDE’s in classical mechanics [26]. 

In this paper, we want to study what is going on with these well-known PDE’s when we 
impose nonholonomic constraints on the rigid body motion. Some partial solutions of these new 
PDE’s are discussed in literature, see, e.g., [14, 20, 29] and references within. Our main aim is 
to prove that these PDE’s for the nonholonomic Chaplygin, Suslov and Veselova systems can be 
easily reduced to the well-studied PDE’s for the Hamiltonian vector field (1.1). Consequently, 
we can directly obtain all the possible integrable perturbations of these nonholonomic systems 
directly from the well-known integrable potentials of the Hamiltonian mechanics. 

The necessary references to the main aspects of nonholonomic mechanics can be found in 
several papers of Sergio Benenti dedicated to the analysis of nonholonomic mechanical systems 
[3, 5, 6 , 7]. Following in the steps of these papers we will consider nonholonomic systems using 
only the knowledge of the basic notions of analytical mechanics, i.e., utilizing a ‘user-friendly’ 
approach to the dynamics of nonholonomic systems proposed by Sergio Benenti [5]. 

This paper is organized as follows. In Section 2 we will introduce Bertarnd-Darboux equation 
for the holonomic particle on the plane and its counterpart for the holonomic particle on the 
sphere. Section 3 contains the main result of integrable perturbations for the Suslov system. 
We will show that integrable potentials for the nonholonomic Suslov problem satisfy to the 
standard Bertarnd-Darboux equation for the holonomic particle in the plane. Sections 4 and 5 
are devoted to the review of the known integrable potentials for the Veselova and Chaplygin 
systems. We will give the Bertrand-Darboux type equations for these systems and will show 
how these equations are reduced to the Bertarnd-Darboux equation for the holonomic particle 
on the sphere. Section 6 deals with two nonholonomic systems on the plane and contains the 
Bertrand-Darboux type equations for the nonholonomic oscillator and the Heisenberg system 
(nonholonomic integrator). Finally, we briefly discuss the quasi-integrable potentials for the 
Suslov and Veselova problems, which were introduced by Llibre et al. 
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2 Integrable potentials on the plane and sphere 


In [9] Bertrand studied the Newton equations for a particle on the plane 


d 2 qi 
dt 2 


= F, 


d 2 q 2 
dt 2 


= F 2 , 


Fu = - 


9V(qi,q 2 ) 

dqk 


( 2 . 1 ) 


and tried to solve equation (1.4) using linear, quadratic or fractional (linear/linear) anzats in 
the velocity for the additional integral of motion. 

In particular, according to Bertrand [9], if 


2 

H 2= KijiQi’QdQiQj + U{qi,q 2 ), 
i,j =1 


then one equation H 2 = 0 yields two systems of PDE’s. The generic solution of the first system 
of equations for the coefficients K t] (q\ , q 2 ) 


H 2 = (-|«2 - P2Q2 + p\ + (-|?i - PiQi + P2 


2 ) 

+ faoi 09 + Bt 09 + Bnai 






depends on the six constants of integration a, (3\, (32, 711 , 712 , 722 - 

In order to describe the forces F\ 2 that should act on the particle Bertrand also extracted 
one linear second-order PDE on potential V from the coupled system of equations on the poten¬ 
tials V and U 


(aqiq 2 + fiiqi + (3 2 q2 + 712 )(^^ “ ^ 11 V) 

+ (aql - aq 2 + 2(3^1 - 2/3 2 <?2 + 7n - 722)<9 i 2 E 

+ 3(aqi + (3i)d 2 V - 3 (aq 2 + ^(diV = 0, (2.2) 

where d{ = d/dqt and dik = d 2 /dqjdqk- 

In [9] Bertand studied only some partial solutions of this equation, whereas in [17] Darboux 
gave a complete solution and, therefore, now equation (2.2) is called the Bertrand-Darboux 
equation. Later on Darboux results were included almost verbatim in the classical text of 
Whittaker on analytical mechanics, see historical details in [32], 

Ideas used by Bertarnd and Darboux to solve the Bertrand-Darboux problem were generali¬ 
zed to study Hamiltonian systems defined in Euclidean spaces of higher dimensions and in other 
(pseudo) Riemannian manifolds. According to Eisenhard [21], the first system of PDE’s is the 
Killing equation for the Killing tensor of second order with vanishing Haantjes torsion, hereafter 
called the characteristic Killing tensor on the Riemannian manifold Q. The second system of 
PDE’s on the various Riemannian manifolds was studied by Kalnins, Miller [25], Benenti [2, 4] 
etc. 

In particular, Benenti formulated the following proposition. 

Proposition 1. A natural Hamiltonian H = ^2gijPiPj + V on the cotangent bundle T*Q of 
a Riemannian manifold Q is separable in orthogonal coordinates iff on Q there exists a Killing K 
of second order with simple eigenvalues and normal eigenvectors, so that 

d(KdV) = 0. (2.3) 

Separable Hamiltonian flow has a necessary number of first integrals, which can be directly 
calculated from the characteristic Killing tensor K, which satisfies to the Killing equation 


V a K]g 7 + V pF-ya + V 7 A ai a — 0, 
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where V is the Levi-Civita connection of the Riemannian metric [ 2 , 4 ], Tensor K has normal 
eigenvectors if and only if its Haantjes torsion is equal to zero. Integrable systems associated with 
the Killing tensors of second order with nontrivial Haantjes torsion were found only recently [ 38 ]. 


2.1 Bertrand—Darboux type equation on the sphere 

Let us consider the standard Hamiltonian vector field describing rotation of a rigid body fixed 
at the point 

M = M x uj, 7 = 7 x uj. ( 2 - 4 ) 


Here M is the angular momentum, uj = AM is the angular velocity, 7 is the constant unit vector 
in a moving frame an A is the diagonal inverse to I matrix 

(a\ 0 0 \ 

A = I -1 = I 0 a 2 0 , a k = ~. (2.5) 

\0 0 a 3 J 4 

According to Euler, there are two first integrals of second order in momenta 
H 1 = ^(M,AM), H 2 = (M,M) 
and two geometric integrals 

Ci= 7i + 72 + 73; c 2 = J 1 M 1 + 7 2 M 2 + 73 M 3 , 


which are the Casimir functions of the underlying Poisson structure. 

Let us consider perturbation of the free motion ( 2 . 4 ) by adding forces associated with potential 
field Vi (7). For the perturbed Hamiltonian vector field 


M = M XW + 7X 


dV\ (7) 
<9 7 ’ 


7 = 7 x uj, 


( 2 . 6 ) 


functions C\ 2 are also constants of motion and, therefore, we can exclude the third component 
of the Poisson vector 7 from other calculations 


73 = yjl -7i -72- 

Substituting the polynomials of second order in momenta 

Hi = ^(M, AM) + V\ ( 7 ), H 2 = (M, M) + V 2 ( 7 ) (2.7) 

in the equations H\ 2 = 0 one gets three partial differential equations on Vh ?2 ( 7 ) 

3 2 (aiK 2 - 2 Vi) = 0 , <9i(a 2 V 2 — 2Vi) = 0 , 

03 ( 72^1 V 2 - 7i<9 2 V 2 ) + 27i<9 2 Vl - 272^1 Vi = 0. (2.8) 

Here Vi ) 2 ( 7 ) are functions on two independent variables 71 , 72 and d/d'jk = dk- 
The generic solution of these equations (2.8) 

Vi = (o 2 a 3 7i + aia 3 7| + 010273), V2 = ai 7 i + a 2 ll + °373 

is associated with the Clebsh system [26]. 
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At C2 = (7, M) = 0 the phase space is equivalent to the cotangent bundle of the two- 
dimensional sphere [ 12 ]. In this case equation H2 = 0 yields only two equations 

(72(^2 - a 3 ) - a 2 )< 9 i V2 - 7172(01 - a 3 )d 2 V 2 + 2 d\V\ = 0 , 

(71(01 — a 3) — 01)^2 V2 — 7172(02 — 0,3) < 9 1 ^2 + 2 S 2 Fi = 0 , ( 2 . 9 ) 

from which we can easily get the equation on one potential. Namely, when equations are 

differentiated by 7^2 and subtracted from each other, the result is 

3i((7 2 2 («2 - 03) - a 2 )diV2 - 7172(01 - a 3 )d 2 V 2 ) 

- ((71(01 - o 3 ) - ai)d2V 2 - 7172(02 - a 3 )d 3 V 2 ) = 0. (2.10) 

It is well-known that the characteristic equation ( 2 . 3 ), ( 2 . 10 ) has a continuum of solutions 
labelled by two arbitrary functions Gip 

w _ 02^1(01) - U\G 2 {U2) _ Gi(ui) - G 2 (u 2 ) 11 1 

*1 r >/ \ 5 *2 • ( 2.11 J 

2(u 2 - Ui) U2 - Ui 

Here u\, u 2 are sphero-conical coordinates on the sphere 


7 < = 


(Ul - Oj)(-U2 
(Oj flj)(o. m 





If p Ul 2 are the corresponding momenta defined by relations 

Mi = 2eij - 7J7m( ° j -° m ) ( ( a, - ai)p» - (a, - u 2 )p„), 

1*1 - 02 

then first integrals Hi ,2 (2-7) satisfy to the following separation relations 

4 (rtj - ai)(u* - a 2 )(ui - a 3 )p 2 u . + Gifa ) - UiH 2 + 2Hi = 0, * = 1,2. (2.12) 

If we take homogeneous polynomials of IV-th order G\ >2 = u N , the corresponding poten¬ 
tials Vi^ satisfy to the well-known recurrence relations 

V} 1] = 0, H 2 (1) = 1, 2V} N) = pVi N ~ 1] , P 2 (7V) = aV 2 {N ~ 1] - 2H 1 (7V_1) , 


where 


cr = (u\ + u 2 ) = (ai + a 2 + a 3 ) - a ±yf - 027 ! - “ 373 ; 

p = uiu 2 = a 2 a 3 7 i + «ia 372 + 010273- 

According to Bogoyavlenskii [12] these potentials are equal to 



(2.13) 


(2.14) 


Here [z] is an integer part of the rational number z. If G\ i2 (u) = u K , one gets rational potentials 


(-K) _ ~ 1 /u£ __ 


,K 


V.\ 


uJ< - uf 


Vn 


(K) 


(U\U 2 ) K U 2 ~ U\ 


rj< 


U 2 ~ U\ 
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and 


(-K) _ u 2 /u f -ui/u$ _ 


,K 


2 V, 


u, 


K+l 


— U 


K+l 

1 


U 2 ~ U\ 


(uiu 2 ) 


K 


U 2 - U\ 


T/ (^+l) 

v 2 


Of course, any linear combination of these polynomial and rational potentials also satisfies (2.9). 
For instance, at N = 2 one gets the Neumann system 


h i 2) = AM ) ~ fl 2«37l - “10372 - ai«273> 

' = (M, M) + ai7i + a 275 + ®373) 


and the Braden system at K = 1 
H 


(-1) = I^ M AM ) + Q !7i + Q272 + q 373 - (7i + 72 + 73 ) ( Q 1 + «2 + as) 
1 2 v 020371 + «1«372 + a l«273 

H { ~ 1] = (M, M) + 1 


02037^ + Oia 3 7| + 01027!' 


In [41] Wojciechowski presented another family of potentials associated with other symmetric 
functions Gy 2 (2.11) on variables u\ t2 and parameters ai, 02 , 03 . For instance, if 


(u) = u(u 2 - (01 + a 2 + a 3 )u + oio 2 + oia 3 + a 2 a 3 ), 
( 3 ) 

F\ 2 (u) = u(u - ai)(u - a 2 )(u - a 3 ) 


then the second integrals of motion read as 

H 2 ] = (M, M) + Y a hi ~ (X ai7 * 2 ) 5 

H 2 ] = (M, M) + Y a hi - 2 (X) a ^f) (X! + (X ar/ f) ■ 

More complicated functions G +2 (2.11) yield more complicated potentials, for example, rational 
functions 


Gi- 2 (u) = Y + 


(oi - a 3 )(ai - a 2 )bf (a 2 - a 3 )(a 2 - ai)bj (a 3 - ai)(a 3 - a 2 )b\ 


Z— 1 


u — a\ 


+ 


u- a 2 


+ 


u- a 3 


give rise to the Rosochatius potentials [31] 


H 


(Ros) _ 1 


1 


= -(M, AM) + 


^1(0273 +« 37 l) . ^i( a l 73 + a 37 l) . 63(017! + 027!) 


7i 

Ff( Ros ) = (M, M) + ^| + t| + t| 


+ 


72 


+ 


7s 


72 


73 


up to the constant terms, see also [24], 

In [19] Dragovic considered functions G 1,2 (2.11), which are the Laurent polynomials, and 
ingeniously coupled the corresponding potentials V\ 2 together with the Appell hypergeometric 
function. 

These familiar well-studied potentials are permanently rediscovered both in holonomic [24, 39] 
and nonholonomic mechanics [22, 29]. 
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3 Suslov problem 


One of the most widely known mechanical nonholonomic systems is the Suslov problem descri¬ 
bing motion of a rigid body under the following constraint on its angular velocity 

(u, a) = 0, (3.1) 


where a is a fixed unit vector in the body frame [33]. It means that there is no twisting around 
this vector a. 

Imposing this constraint we have to add some terms with the Lagrangian multiplier to the 
initial Hamiltonian vector field 


IcL? = Ixcu + 7X 


dVi{i) 


+ A a, 


7 = 7 x u. 


(3.2) 


Differentiating the constraint (3.1) by time and using the equation of motion we obtain 


A 


1 

(l _1 a,a) 



1 

X 0 J + -7 X 


dVihY] 
<97 J- 


Vector field (3.2) preserves the mechanical energy 


H\ = (M,AM) + Vi(rf) 


(3.3) 


and the geometric constants of motion 


Ci = ( 7 , 7 ) = 1 , C 2 = (w, a) = 0, 


which allows us to remove the redundant variable from the calculations 
73 = ^1-7i-72- 

If we assume that a is an eigenvector of the tensor of inertia [14, 20, 29], i.e., that tensor of 
inertia is diagonal 

fh 0 0 \ 

1=0 h 0 , 

\0 0 h) 

and vector a is equal to a = (0,0,1) in some coordinate frame, the constraint is trivial 


w 3 = 0. 


Substituting the standard anzats for the second integral of motion 

H 2 = /1 (7)^1 + / 2 ( 7 )w| + h{l)^im + V 2 {i) 

with unknown functions f). and V 2 on the components 712 of the Poisson vector in the equation 
H 2 = 0 we obtain the well-known Bertrand-Darboux equation on potential [9, 17]. 

Proposition 2. For vector field (3.2) the following statements are equivalent: 

1. There is an additional independent first integral of second order in velocities 

f a\fhh 2 . o rr , 7nv / 7A 2 

Hi = ( 2 ^1 + W'm + ) “1 

f ay/hh 2 , 0 / 7 “ , l 22 \/h\ 2 

+ [-^—72 + Mhl2 + -^=-J w 2 

+ [ay/I1I2T1T2 + /3iV^i72 + foy/hli +712) <^iu;2 + ^2(7)- 
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2. Potential V\ satisfies the Bertrand-Darboux equation (2.2) with 

Qi = q 2 = 72 (3.4) 

3. Potential V\ is separable. A characteristic coordinate system for the Bertrand-Darboux 
equation provides separation for V\ and can be taken as one of the following four orthogonal 
coordinate systems on the qi^-plane: elliptic, parabolic, polar or Cartesian. 

The proof is completely similar to the one for the original Bertrand-Darboux theorem [17, 32]. 

This result allows us to suppose that the nonholonomic Suslov system is equivalent to the 
holonomic motion on the plane with coordinates q \, q 2 after some singular change of time, but 
its study is out of the framework of the present note. 

Another reduction to the Bertrand-Darboux equation was proposed in [29]. Integrable vector 
field remains integrable for any fixed value of mechanical energy Hi, for instance on the zero- 
energy hypersurface 

H\ = 0 . 


In the Hamiltonian case the separation of variables of this null Hamilton-Jacobi equation is 
equivalent to the ordinary separation of the image of the original Hamiltonian under a generalized 
Jacobi-Maupertuis transformation [ 8 ]. 

If we substitute potential 


V ' = -2 


and velocities 
mu 2 


mui 


h 


LOl = 


h 


+ 


mun 


CO 2 —- — 


fil 

h' 


C03 = 0 


from Theorem 1 in [29] into the first integral H\ (3.3), one gets Hi = 0. Here p,i :2 are functions 
on the components of the Poisson vector 71 2 , which satisfy equations (3.2) 


/1 2 2^1 • /1 2 2^2 

71 = -y 1 - 7i - 72 72 = -y 1 - 7i - 72^- 

Differentiating these equations with respect to time we obtain the Newton equations (2.1) 
on qi t 2 (3.4) with forces labelled by two functions Hi )2 {qi,q 2 ). Thus, on the zero-energy hy¬ 
persurface of the initial Hamiltonian one gets an initial Bertrand problem with the well-known 
solution. Of course, on this zero-energy hypersurface we can find other solutions associated 
with nonitegrable potentials on the whole phase space. In [29] such potentials were called quasi- 
implicitly integrable or locally integrable potentials. 


4 Veselova system 

Let us consider the nonholonomic Veselova system describing the motion of a rigid body under 
the following constraint 

(w, 7 ) = 0, (4.1) 

where 7 is a unit Poisson vector fixed in space [40]. It means that there is no twisting around 
vector 7 . 
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According to [40] this constraint shifts the initial Hamiltonian vector field (2.4) 

M = M x io + Ay, 7 = 7 X 141 , (4-2) 

where the Lagrangian multiplier A is chosen so that the constraint (4.1) is satisfied at any time 

(AM x M, A 7 ) 

(Ay, 7) 

There are integrals of motion of second order in momenta 
Hi = ^(M, AM), H 2 = (M, M) - (y,y)- 1 (y,M ) 2 
and two geometric constants of motion 


Ci = (7,7), C 2 = (y,u;) = 0. 


In the presence of the potential field equations of motion (4.2) become 

dW\ ( 7 ) 

M = Afxw + Ay + 7X —--, 7 = 7 x 0 ;, 


where 


A = 


<9y 


(AM x M + 7 x dWi( 7 )/ 07 , Ay) 


(4.3) 


(A 7 , 7 ) 

As usual, functions C\p remain constants of motion, and we can exclude the redundant variable 


73 = y 1 - 7? - 71- 

Vector field (4.3) is a conformally Hamiltonian field, see, for instance, [15, 36]. Substituting the 
following anzats for integrals of motion 

H 1 = *(M, AM) + HA ( 7 ), H -2 = (M, M) - ( 7 , y)- 1 (y, M ) 2 + W 2 (y) (4.4) 


in H\ )2 = 0 one gets two first-order equations on potentials W\. 2 

2(72 (°2 1 — a 3 1 ) _ a 2 l )diW\ - 27172(07 1 - af l )d 2 W\ + d\ W 2 = 0, 

2(71 (a ] -1 - aj x ) - - 27172 ( 0^ 1 - l )d{W\ + d 2 W- 2 = 0. (4.5) 

Proposition 3. After the inversion of parameters ak —X 07 1 an< ^ substitution 
IP 2 = 2Vl, Wi = y 

equations (4.5) coincide with equations (2.9) for potentials on the two-dimensional sphere. 

Thus, all the integrable potentials for the nonholomic Veselova system are easily expressed via 
well-known integrable potentials Vi 2 f° r the holonomic system on the two-dimensional sphere. 
For functions G\ 2 (2.11) which are the Laurent polynomials in u\ i2 these expressions were found 
in [20]. We only want to note that it is true for any integrable potentials. 

Following Theorem 2 in [29] let us fix the values of velocities by equations 

CLW 172 - / 2 W 271 - T 2 = 0, poo 3 — A>i = 0, uqyi + W 272 + W 373 = 0, 
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where p = Jhhh (tt + 77 + S) 
and potential 


and ^ 1 , 2 ( 7 ) are functions on 7 . Substituting these velocities 


Wi = 


+ 'I'l 

2(/i72 + J 27 1 ) 


into the mechanical energy Hi (4.4) one gets H\ = 0. 

The remaining three equations of motion of the components of the Poisson vector 7 are easily 
reduced to equations of motion for the holonomic particle on the sphere with forces labelled by 
two functions ^ 1 , 2 ( 7 )- Thus, on the zero-energy hypersurface of the initial Hamiltonian one gets 
standard characteristic equation ( 2 . 10 ) with the well-known solutions. 


5 Chaplygin ball 


As in [16] we consider the rolling of a dynamically balanced ball on a horizontal absolutely rough 
table without slipping or sliding. ‘Dynamically balanced’ means that the geometric center 
coincides with the center of mass, but mass distribution is not assumed to be homogeneous. 
Because of the roughness of the table this ball cannot slip, but it can turn about the vertical 
axis without violating the constraints. 

After reduction [16] motion of the Chaplygin ball is defined by the following vector field 
M = M x lo, 7 = 7 x oj. (5.1) 


Here M is the angular momentum of the ball with respect to the contact point, c o is the angular 
velocity vector of the rolling ball. Its mass, inertia tensor and radius will be denoted by rn. 
I = diag(/i, I 2 , 13 ) and b respectively. All the vectors are expressed in the so-called body frame, 
which is firmly attached to the ball, and its axes coincide with the principal inertia axes of the 
ball. 

The angular velocity vector is equal to w = A g M, here matrix 
A g = A + dg(y)A 7 <g) 7 A 


is defined by the nondegenerate matrix A (2.5) and function 

1 


g(7) = 


d = mb 2 . 


(5.2) 


1 - d { a^l + 027I + CI37!) ’ 

It is easy to prove that vector field (5.1) preserves two polynomials of second order in momenta 

H\ = ^ (M, A g M ), H 2 = (M, M) 

and two geometric constants of motion 

Ci = 71 + 7! + 73 = 1 , C 2 = 71 Afi + 72^2 + 73M3, 


see details in [13, 34, 37]. 

Indeed, equations of motion of the ball in the potential field 


M = M XCJ + 7 X 


^1(7) 

97 : 


7 = 7 x 10 


(5.3) 


have the same form as the equations (2.6) in rigid body dynamics. In fact, the principal difference 
between holonomic and nonholonomic systems is hidden within the relation of the angular 
velocity to the angular momentum. 
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According to [35, 36], integrals of motion for the Veselova system are expressed via the 
integrals of motion for the Chaplygin ball. So, we can easily express integrable potentials for 
the Chaplygin system via integrable potentials for the Veselova system W\ 2 (4.5) and then via 
integrable potentials Vi ,2 (2.9) for the holonomic system on the two-dimensional sphere. Of 
course, we can get the same result by directly substituting standard ansatz 

= i(M, A g M) + Vi( 7 ), H 2 = (M, M) + C 2 ( 7 ) (5.4) 

in Hi )2 = 0 we obtain 

ai{a 2 - a 3 ) 7 i 7 20 iV 2 + ai((a 2 - a 3 ) 7 | + a 3 - 1 )d 2 U 2 + 2g _ 1 0 2 Ci = 0, 
a 2 ((ai - a 3 ) 7 f + a 3 - l)0i U 2 + a 2 (ai - a 3 ) 7 i^ 2 d 2 U 2 - 2g~ l d 3 Ui = 0, 

0372((01 - a 2 ) 7 i + a 2 - l)diU 2 + a 37 i((ai - a 2 ) 7 | - ai + l) 0 2 V 2 
+ 2 g _1 ( 7 i 0 2 Vi - 72 <9iVi) = 0 . 

This system of equations has only one solution 

Hi = "(020371 + “io 3 72 + 01027 !), U 2 = ai'yl + a 2 ^l + a 37 |, 

which coincides with the single solution of the initial system (2.8) associated with the Clebsch 
model. This integrable potential has been found in [27]. 

At C 2 = 0 conditions H\ 2 = 0 are thus 

g 7 i 7 2 (a 2 (ai - a 3 ) 7 ! + ai(a 2 - a 3 ) 7 | - aia 2 + (01 + o 2 - l)a 3 )diU 2 

- g(ai(xf + x\ - 1) ( a 2 xl - l) - a 3 (ai(x\ - 1) + x\(a 2 xl - l)))d 2 U 2 

- 2 7 i 72 0 iVi - 2 ( 7 | - 1 ) 0 2 Vi = 0 , 

g(a 2 (xf + x 2 - l) (a\x\ - 1 ) — a 3 (a 2 (xf - l) +x$(a 1 xj - l))) 0 i U 2 

- g 7 i 72 ( 02(01 - a 3 ) 7 ! + ai(a 2 - a 3 ) 7 | - aia 2 + (ai + a 2 - 1 )a 3 )d 2 U 2 

+ 2(xf - l)0iVi + 2 7 i 7 202 Vi = 0. (5.5) 

Here g = g( 7 ) is the function defined by (5.2). If we change the parameters 

ai a 2 a 3 

ei — 1 -> e 2 — --, e 3 — -- 

1 — 01 1 — ci 2 1 — a 3 

and substitute in (5.5) 

2 Hi = (e 2 e 37 ! + eie 37 | + eie 2 7 3 )C 2 + 2 Vi, 

U 2 = d( 1 + (ei + e 2 + e 3 ) - e\jf - e 2 ^\ - e 37 f + e 2 e 37 ! + eie 37 | + eie 2 7 |) V 2 , (5.6) 

then the equations (5.5) become 

( 7 l ( e 2 - e 3 ) - e 2 )0iVi - 7172 (ei - e 3 )0 2 V 2 + 20iVi = 0, 

( 7 i (ei - e 3 ) - ei)0 2 Vi - 7172 (e 2 - e 3 )0iV 2 + 20 2 Vi = 0. 

It is easy to see that this system coincides with the initial system of equations (2.9) defining 
integrable potentials on the sphere up to oa, —> e^. 

Thus, for the Chaplygin ball imposition of the nonholonomic constraint leads to deformation 
of potentials (5.6) and to replacement of parameters —> e^. 
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Proposition 4. At C 2 = 0 conformally Hamiltonian vector field (5.3) has two integrals of 
motion (5.4) with potentials 

2U\ = pV '2 + ‘IV \, U 2 = d(p + a + 1)V 2 . 

Here V\p are integrable potentials on the sphere (2.11) after replacement of parameters a & —> ek, 
and 


o = (7i + 72 + 73) ( e i + e 2 + e 3 ) - ei7? + e 2 y| + e 37 3 , 

P = e 2 e 3 7i + e i e 37 2 + eie 2 7 3 

are £/ie same polynomials of second order in variables 7 as above (2.13). 

Following to S.A. Chaplygin [16] we can introduce the sphero-conical coordinates u\, u 2 


H = 


(u\ - ej)(u2 
(&j Cj)(e m 



* + j + m, 


and explicitly present some solutions of the Bertrand-Darboux type equations (5.5). Namely, 
integrals of motion H±p (5.4) satisfy the separation relations 


4(ei - Ui)(e 2 - Ui)(e 3 - uf) N 

d(e 1 + l)(e 2 + l)(e 3 + l) Pui + ,( “ iJ 


UjH-2 

d(ui + 1 ) 


+ 2H X = 0, 


* = 1 , 2 , 


which can be considered as gentle deformation of the initial relations (2.12). Thus, separable 
potentials in this case read as 


U 2 


Ui 


dfia 1 + 1 )(u 2 + l)(G , i(ui) - G 2 {u 2 )) 
u 2 - U\ 

uiu 2 {Gi(ui) - G 2 (u 2 )) + u 2 Gi{ui) - uiG 2 {u 2 ) 
2 (u 2 - ui) 


The passage to limit d —> 0 reduces equations of motion for the Chaplygin ball (5.3) to the 
standard Euler-Poisson equations. However, at d —> 0 we have to simultaneously change the 
definition of the second potential U 2 in (5.6) and, therefore, we present another family of solu¬ 
tions for equations (5.5). 

Let us introduce variables v\p 


7i = 


(1 - daj)( 1 - da m ) / (vi ~ aj){v 2 ~ a>i) 
(1 - dv i)(l - dv 2 ) V (a,j - ai)(a m - af) ’ 


i + j + m, 


and the conjugated momenta p tJl 2 , see [34] for details. In this variables the separated relations 
have the following form 

4(1 - dvf)(vi - ai)(vi - a 2 )(vi - a 3 )pl. + Ufivi) + ViH 2 - 2 Hi =0, i = 1,2, 


and integrable potentials 


U 2 


GM - G 2 (v2) 

v 2 — V\ 


_ v 2 Gi(vi) - ViG 2 (v 2 ) 
' 2 (u 2 - vi) 


are the same functions on variables ui j2 as the integrable potentials on the sphere ( 2 . 11 ). 
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Proposition 5. At C -2 = 0 vector field ( 5 . 3 ) has integrals of motion of second order in veloci¬ 
ties ( 5 . 4 ) if potentials U\p have the same form as integrable potentials on the sphere ( 2 . 14 ) 



Of course, any linear combination of these polynomial and rational potentials also satisfies the 
equations (5.5). These potentials differ from (2.14) by replacement of polynomials a and p for 
the following functions 

? = g(7)(^ + d (a i(a 2 + 03 ) 7 ? + a 2 (ai + 03 ) 7 ! + 03(01 + 02 ) 73 ))* 

Q = §(7)^ + ^010203(71 +72 +73)) > 

which at d = 0 become initial polynomials (2.13). 

For instance, at N = 2 we have the following analogue of the Neumann system 

H[ 2) = i (M, AM) + g, H (2) = (M, M) + ? , 

and at K = 1 the following counterpart of the Braden system 

H[~ 1] = i H^~ l) = (M, M) - -. 

A 6 g 

Of course, we can also single out other families of solutions of the equations (5.5), for instance, 
see [20]. 

6 Nonholonomic oscillator and Heisenberg system 

Let us consider the Lagrangian of the particle in Euclidean space M 3 

L = T j{x 2 + y 2 + z 2 )- v i x ,y,z), (6.1) 

where m is the mass of the particle. When this system is subject to a nonholonomic constraint, 
the resulting mechanical system may or may not preserve energy and the phase space volume, 
and their integrability and reduction theories are completely different from the Hamiltonian 
case [1, 18, 23]. In this Section we consider two first-order nonholonomic constraints, which 
displays all the basic properties of first-order nonholonomic systems in the control theory [11]. 

The first nonholomic constraint and potential in (6.1) for the so-called nonholonomic oscil¬ 
lator have the following form 



f = z - kyx = 0, 


k eR 
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whereas second constraint and potential in (6.1) for the so-called Heisenberg system read as 


f = z — (yx — xy) = 0, V = (1. 


The Heisenberg system (nonholonomic integrator) can be pointed out as a benchmark example 
of nonholonomic system with a first-order nonintegrable constraint, which mimics the kinematic 
model of a wheeled mobile robot of the unicycle type. 

In generic case at V ( x, y, z) = V(x, y) in (6.1) the third degree of freedom also decouples from 
the rest of the system and after nonholonomic reduction we obtain a two-degrees of freedom 
system of the Chaplygin type [11, 23, 30]. For the such generalized nonholonomic oscillator the 
reduced equations of motion are 

i= m’ Px = ~ 1 + fcV i^yPx'Py + d x V), p y = ~d y V. (6.2) 

For the generalized Heisenberg system the reduced equations of motion read as 


Px 


Py 


x = —, y = 
m m 

„2 


Px = 


0 x 2 + l)d x V + xydyV 
m{ 1 + x 2 + y 2 ) 


Py = 


(: y 2 + 1 )dyV + xyd x V 
m( 1 + x 2 + y 2 ) 


(6.3) 


Below we will study potentials V(x,y) in (6.2) and (6.3), so that the corresponding four¬ 
dimensional vector fields X have an additional first integral and possess an invariant volume 
form. Of course, equations on these potentials have the form of the characteristic equation (2.3) 
and can be considered as an analogue of the Bertrand-Darboux equation (2.2). 


6.1 The generalised nonholonomic oscillator 

The vector field for the reduced nonholonomic oscillator (6.2) after the following change of 
variables 


Pi 


= \/k 


2 y 2 + 1 p x 


P2 = 


Py 


\/k 2 y 2 + 1 


q i = x, 


<?2 = y 


becomes the conformally Hamiltonian vector field 


X = -pPdH 1 , P =(° j o) 


with respect to the canonical Poisson bivector P and reduced Hamiltonian 
2 


Hi = ^2 SijPiPj + V(qi,q 2 ) 
i,j=l 


P \ , pl(k\l + 1 ) , , 

2 ^ + -2m- + 


Conformal factor 

1 

M V k2 qi + 1 

is a nowhere vanishing smooth function on an open dense subset of the plane q 2 7 ^ 00 , which 
defines an invariant volume form H = ydq A dp. 

Substituting a linear function in velocities 


H 2 = gi(qi,q 2 )pi + 52 ( 91 , <?2 )P 2 
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into the equation iP 2 = 0 one gets 

gi = ci + c 2 In (kq 2 + \Jk 2 ql + l), g 2 = - \Jk 2 ql + l(c 2 fcgi - c 3 ) 


and 


ci ln(/cg 2 + y/k 2 q$ + !) ( kqj fln(kq 2 + y/k 2 q 2 + l)\ \ 

V ^ = G [ -*- C2 ^ + i—y^i + i 

If we want to consider a single valued integral H 2 , we have to put C 2 = 0 and 

c 3 = 0, V = G(qi) or a = 0, V = G(q 2 ). 

Substituting polynomials of second order in velocities 
2 

H 2= K ij(^^ 2 )PiPj + U(qi, q 2 ), 

0=1 

where Kjj and U are single valued functions on gi, q 2 in the equation H 2 = 0, we obtain the 
following expression for the second integral of motion 

ci {k 2 ql + kq 2 y/k 2 q\ + 1 + l) 2 2 , n 2 , TTl ^ 

H 2 = -;- 7 = 9 -P 1 P 2 + c 2 {k q 2 + 1 )p 2 + U{qi,q 2 ), 

kq2 + \A ^ + 1 

and the following counterpart of the Bertrand-Darboux equation 


ci 


■\{k 2 ql + l)d 2 2 V + k 2 q 2 d 2 V - duV) - 2c 2 9i 2 P = 0. 


y/Wq^+l 

This equation has one physical and one formal solution 

ci = 0 , V = G 1 (q 1 ) + G 2 (g 2 ) and c 2 = 0 , V = G 1 (q+) + G 2 (g_), 

where 

ln(fcg 2 + ^/k 2 qf+l) 


(6.4) 


q± = 91 ± 


jfe 


Thus, for the nonholonomic oscillator we obtain only trivial perturbations in the framework of 
the Bertrand-Darboux method. 

6.2 The generalized Heisenberg system 

The vector field for the reduced Heisenberg system (6.3) after the following change of variables 


Pi = 


m((l+y 2 )p x - xypy) 


1 + x 2 + y 2 
is conformally Hamiltonian vector field 

X = -/rPdHi, P= ' ° 1 


P 2 = 


m((l + x 2 )p y - xyp x ) 
1 + x 2 + y 2 


Q i = x, q 2 = V 


V -I 0, 

with respect to canonical Poisson bivector P and reduced Hamiltonian 

2 ' „2 _|_ '2 _|_ j 

Hi SijPiPj + V (gi, 92) = 1 2r ^ -(pi + pi + (91P1 + 92P2) 2 ) + P(gi, g 2 ). 

*d=i 


(6.5) 


2m 
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Conformal factor 

/*=(i +<R+<&r l 

is a nowhere vanishing smooth function on an open dense subset of the plane q\ 2 7^ 00 , which 
defines an invariant volume form = pdq A dp. 

Substituting linear function in velocities 

#2 = 01 ( 11,12 )pi + 52 ( 11 , 12)02 

into the equation H -2 = 0 one gets the following first integral 

H 2 = (piqf + 021102 +Pl)ci + (piqiq2 +P2I2 +P2)C2 + (pil2 ~P 2 qi)c 3 , c k G M, 

and potential 

v = G ( (cm + C 2 I 2) 2 + 2c 3 (ciq 2 ~ C 2 I 1 ) + cj + cj - c| \ 

V ( c ll2 - C 2 I 1 - C3) 2 ) 

depending on the arbitrary function G. 

Substituting polynomials of second order in velocities 

2 

H 2= ^2 12 )PiPj + U(q 1 ,q 2 ) 

i,j =1 

into the equation H 2 = 0 one gets the following expression for the second integral of motion 

H2 = (P1I2 -P 2 li)(pi 1 i +P2I1I2 +Pi)ci + (piqj +P2I1I2 +Pi)(pi 1 il 2 + 02 li +P2 )c 2 
+ (Pll2 -P2ll)(pi1ll2 +P212 +P2 )c 3 - (PH2 -P2ll) 2 C4 
+ (l!(li + l)pi + 2(/2li P1JJ2 + (12 + 1 ? + 2 g| + l)po ) c 5 
+ ((li + 2 gi + 12 + l)p 2 + 211I2P1P2 + 1 i (l| + l )^) c 6 + t^(ii, i 2 ). 

In this case equation (2.3) looks like 

Ad u V + 2Bd u V + Cd 2 2V + -- \ -o ( ad i V + bd 2 v ) = 0, 

I + I 1 +I 2 

where A, B, C are the polynomials of second order in q\p 

A = (q'f + l)(iici - c 2 ) + 12 ( 1 ? - l)c 3 + 2gig 2 (c 6 - c 4 ), 

B = 12 (li + l)ci + gi (12 + l)c 3 + (li — 1 2 ) c 4 — (li + l)c5 + (i| + l)c6, 
c = (if + l)(l 2 C 3 + c 2 ) + 11 (li - l)ci + 2gig 2 (c 4 - c 5 ), 

and a, b are the polynomials of fourth order 

a = (2 qf + 2g^l2 + 5gi — I 2 + 3)ci + 11 (li — 3gf + l)c 2 + 2gi g 2 (if + if + 3)c 3 

- 212 ( 1 ? + i| + 3 )c 4 + 4g 2 (li + 1)C5 - 212 (li - if - l)ce, 
b = 2 ^ 152(11 + I 2 + 3) ci + g 2 (3i? — i| — l)c 2 + (2 giif + 2g| — li + 5gf + 3)c 3 

+ 2gi(g? + i| + 3 )c 4 - 2gi(gi - gf + l)c 5 - 4gi(g| + l)c 6 . 

Following Darboux [17] we can find the canonical form of the corresponding Killing tensor 
and a few families of solutions to this equation. For instance, if cq = 1 and other constants of 
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integration are equal to zero, then solutions of the equation H 2 = 0 are labelled by two arbitrary 
functions G\p 


V(qi,q 2 ) = 


gi + 92 + 1 

2m 


G lfe) + - G2 


+ 1 

ql 


U(qi,q 2 ) = qfGi(q 2 ) + G 2 


ql + 1 
ql 


If C 4 = 1 and other constants of integration are equal to zero, then we have solution 


V{qi,q 2 ) = G'i(r) - 


r 2 + 1 


2 mr 2 

associated with polar coordinates on the plane 


G 2 {q>), U(qi,q 2 ) = G 2 (cp) 


qi=rcosp, q 2 = r sirup. 


Here GY 2 are arbitrary functions. In similar manner we can get solutions associated with 
parabolic and elliptic coordinates on the plane, but they are bulky and, therefore, we do not 
present these solution explicitly. 

In order to get these solutions we can also use the Birkhoff method. Namely, let us consider 
a general natural system of two degrees of freedom described in certain generalized coordinates 
by the following Lagrangian 

hj= 1 


According to Birkhoff [10] using change of time t —> r and coordinates (qi,q 2 ) —>• (x,y) this 
Lagrangian can always be reduced to the form of 


L = 



U (x, y). 


Thus, taking the well-known solutions U(x,y) of the classical Bertrand-Darboux equation (2.2) 
and applying the inverse Birkhoff transformation we are able to obtain integrable potentials 
V(q \, q 2 ) for the two-dimensional holonomic system with nonstandard metric. 

We can apply this method for the given nonholonomic case because the corresponding vec¬ 
tor field X (6.5) is a conformally Hamiltonian vector field [30, 36], i.e., it can be reduced to 
a Hamiltonian vector field by changing of time. Recall, that equation H 2 = 0 (1.4) is invariant 
with respect to change of time, which we have to use both in the Birkhoff method and in the 
reduction of the conformally Hamiltonian vector field to the Hamiltonian one. 


7 Conclusion 

In this paper, we consider perturbations of the five well-known two-dimensional nonholonomic 
systems, which are integrable by the Euler-Jacobi theorem. We show that the Bertrand- 
Darboux method is applicable to these systems and that all the obtained Bertrand-Darboux 
type equations in nonholonomic case can be reduced to the Bertrand-Darboux type equations 
in holonomic case. Consequently, we can directly obtain all the possible integrable potentials for 
these nonholonomic systems directly from the well-known integrable potentials of the Hamilto¬ 
nian mechanics. 
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